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ABSTRACT 

The recent development of unified theories of active galactic nuclei (AGN) has in- 
dicated that there are two physically distinct classes of these objects - radio-loud and 
radio-quiet. The primary observational distinctions between the two types are: (1) The 
radio-loud objects produce large scale radio jets and lobes, with the kinetic power of the 
jets being a significant fraction of the total bolometric luminosity. On the other hand, the 
weak radio ejecta of the radio-quiet objects are energetically insignificant. (2) The radio- 
loud objects are associated with elliptical galaxies which have undergone recent mergers, 
while the radio-quiets prefer spiral hosts. (3) The space density of the radio-louds at a 
given optical luminosity is ~ 10 times lower than that of the radio-quiets. Despite these 
differences, the (probably) thermal emissions from the AGN (continua and lines from X- 
ray to infrared wavelengths) are quite similar in the two classes of object. We argue that 
this last result suggests that the black hole masses and mass accretion rates in the two 
classes are not greatly different, and that the difference between the classes is associated 
with the spin of the black hole. 

We assume that the normal process of accretion through a disk does not lead to rapidly 
spinning holes, and propose instead that galaxies (e.g. spirals) which have not suffered a 
recent major merger event contain non-rotating or only slowly rotating black holes. When 
two such galaxies merge, the two black holes are known to form a binary and we assume 
that they eventually coalesce. In the small fraction of mergers in which the two "parent" 
galaxies contain very massive holes of roughly equal mass, a rapidly spinning, very massive 
hole results. It is proposed that such mergers are the progenitors of powerful radio sources, 
in which the radio jets are powered by the spin energy of the merged hole. We calculate the 
distributions of mass and spin for the merged holes from the parent hole mass distribution, 
which is derived from the optical luminosity function of radio-quiet AGN adopting different 
activity patterns. The ratio of the number of radio-loud to radio-quiet AGN's at a given 
thermal (e.g. optical) luminosity is determined by the galaxy merger rate. The required 
fraction of galaxies which merge during the average lifetime (r^ 10^ yrs) of a radio-loud 
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AGN is found to be 10~^, i.e. a merger rate of 1 in ~ 10^ yrs. The Blandford-Znajek 
formalism is then used to predict the radio luminosity and radio luminosity function of the 
merged population. Comparisons between the predicted and observed radio luminosity 
functions constrain the efficiencies with which jet power is extracted from the spinning 
hole and radio emission is produced by the jet. The cosmological evolution of the radio 
properties of the radio-loud objects is related to the increased frequency of merger events 
at earlier epochs. 

Subject headings: black hole physics - galaxies: active - galaxies: interactions - galaxies: 
jets - galaxies: nuclei - galaxies: quasars: general 
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1. INTRODUCTION 

Over the past decade, a considerable unification in our understanding of the phe- 
nomenology of active galactic nuclei (AGN) has been achieved. Because radiation is not 
emitted isotropically by these nuclei, intrinsically similar types of AGN may appear quite 
different when viewed from different directions. Two kinds of anisotropy are recognised. 
One kind is emission in a narrow cone, which is almost certainly a result of beaming by 
a relativistic jet. The second kind of anisotropy involves emission in a wider cone, which 
is a consequence of either shadowing of an intrinsically isotropic, optical-ultraviolet source 
by a thick (or warped) dusty torus or an intrinsically anisotropic optical-ultraviolet emit- 
ter. Radio-quiet objects form one class of active galaxies, with Scyfcrt I's and radio-quiet 
quasars viewed close to "pole-on" (so the broad line region is seen directly) and Seyfert 
2's viewed more equatorially (so the broad line region is hidden by the dusty torus). The 
second class comprises low power (Fanaroff-Riley class I, FRI) radio galaxies and BL Lac 
objects. A BL Lac is seen when the nucleus is viewed within the radiation beam of a rela- 
tivistic jet which fuels the large-scale radio lobes. At larger viewing angles, the large-scale 
lobes dominate the radio emission and the object is a classical FRI radio galaxy. Lastly, 
OVV quasars, radio-loud quasars, broad- line and narrow- line radio galaxies (FRII) appear 
to form a third class, with the sequence representing increasing angular difference between 
our line of sight and the axis of the source. In this class, both kinds of anisotropy are 
believed to be present. These issues have been recently reviewed by Antonucci (1993). 

The reason for the large range in radio powers is unknown. There has, however, been 
considerable discussion concerning whether the radio-quiet and radio-loud objects form 
distinct and separate classes, or whether there is a continuum of radio properties in just one 
class of object. Several arguments favor the former hypothesis, as we now summarise, (i) 
There is a difference in the nature of the host galaxies, at least at low redshift. Radio-loud 
AGN are never found in spiral hosts, but almost always reside in elliptical or elliptical-like 
galaxies (a few are associated with SO's). Conversely, relatively few radio-quiet AGN are 
found in ellipticals, (ii) Imaging surveys have shown that the mean absolute magnitudes 
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of the underlying galaxies of radio galaxies and radio-loud quasars are the same (<My > 
= -23.3; Ho = 50 km s~^ Mpc~^ and qo = are used throughout the paper), consistent 
with the notion that they are identical objects and that the observability of the quasar 
nucleus is determined by variability or aspect dependent effects (e.g. Malkan 1984; Smith 
et al. 1986; Veron-Cetty & Woltjer 1990; see Hutchings 1987 for a contradictory result). 
In contrast, the host galaxies of radio-quiet quasars are 0.6 - 1.0 mag less luminous (e.g. 
Smith et al. 1986; Veron-Cetty & Woltjer 1990; see also review by Woltjer 1990). It has 
further been argued (e.g. Smith et al. 1986; Hutchings, Janson & Neff 1989) that the host 
galaxies of radio-loud quasars are ellipticals and those of radio-quiet quasars are spirals, 
consistent with the distinction at low redshift. These results show that the radio- loud 
and radio-quiet classes inhabit different types of host galaxy, (iii) The radio luminosity 
functions of Seyfert galaxies and radio galaxies are distinct and do not join onto each other 
(Meurs & Wilson 1984). Further, radio- loud AGN comprise an intrinsically rarer class of 
object than radio-quiets (e.g. Padovani 1993). (iv) Radio surveys of optically selected 
quasars commonly show a bimodal distribution of flux densities (e.g. Strittmatter et al. 
1980; Kellermann et al. 1989). Peacock, Miller & Longair (1986) favored the existence 
of two quasar populations with widely differing radio luminosities and differing optical 
luminosity functions. This conclusion was based on their finding that the radio properties 
of optically faint quasars conflicted with a single population model (Schmidt 1970) in which 
the optical and radio luminosities are correlated. In subsequent work. Miller, Peacock & 
Mead (1990) performed a radio survey of optically selected quasars covering a wide range 
of optical luminosities, but only a narrow range of redshifts. These authors found strong 
evidence for the existence of two populations of quasar - one having radio luminosities 
LbGHz > 10^^ WHz~^ and the other having L^qhz < 10^^ WHz~^ , with no quasars of 
intermediate luminosity. They further concluded that the fraction of quasars which are 
radio-loud is a function of redshift and probably optical luminosity also, and that these 
effects may be understood in terms of the mixing of two quasar populations with differing 
luminosity functions and differing amounts of cosmological evolution, (v) In both radio- 
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quiet (de Bruyn & Wilson 1978) and radio-loud (Baum & Heckman 1989; Rawlings & 
Saunders 1991) objects, there is a strong correlation between the steep-spectrum radio 
and the total [OIII] A5007 luminosities These correlations follow roughly parallel loci in 
the radio - [OIII] plane, but with the radio-loud objects a factor of ~ 10^ more luminous 
at 5 GHz than the radio-quiets for a given [OIII] luminosity. Few objects fall between the 
two classes (Rawlings 1993), suggesting a genuine dearth of "radio-intermediate" objects. 

Observations indicate that the properties of radio-loud and radio-quiet objects are 
similar at infrared, optical, ultraviolet and X-ray wavelengths (e.g. Steidel & Sargent 1991; 
Sanders et al. 1989; Neugebauer et al. 1986). The shapes of the local optical luminosity 
functions of radio-loud and radio-quiet quasars are not significantly different (Padovani 
1993), although the former function is subject to large observational uncertainties. There 
is now considerable evidence that most of the infrared - X-ray continuum is dominated by 
thermal emission which, in the context of black hole accretion models, is related to the 
black hole mass and the ratio of the accretion rate to the Eddington accretion rate (i.e. the 
Eddington ratio). These similarities of the thermal properties of the two classes suggest 
that the distributions of black hole masses and Eddington ratios are not widely different. 
It is then reasonable to infer that the distinction between the radio-loud and radio-quiet 
classes is related to the spin of the black hole. 

The hypothesis (Toomre & Toomre 1972) that many elliptical galaxies are formed by 
the merging of spirals is supported by recent observations (Schweitzer 1982; Kormendy 
& Djorgovski 1989). Images of radio-loud AGN show clear evidence for morphological 
peculiarities (e.g. tails, bridges, shells, dust lanes, second nuclei) indicative of a galaxy 
interaction or merger involving at least one disk galaxy (e.g. Heckman et al. 1986; Smith 
et al. 1986; Hutchings 1987). Such evidence for interactions/mergers is found much less 
frequently in radio-quiet objects (e.g. Hutchings & Neff 1990). Thus, there is clear evidence 
that the radio- loud AGN phenomenon is triggered by a recent (< 10^ yrs ago) interaction 
or merger. 

The present paper presents a model which attempts to account for the dichotomy 
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between the radio-loud and radio-quiet classes. We adopt the standard supermassive black 
hole (SMBH) model of nuclear activity (e.g. Rees 1984). In essence, the thermal emission 
(optical-ultraviolet continuum, at least some of the X-rays, and infrared reradiation from 
dust grains) is powered by accretion. On the other hand, the mechanical output of the 
AGN (i.e. the radio jets) is considered to be powered by the spin of the black hole (e.g. 
Blandford 1990 and references therein). We shall, however, assume that the normal process 
of accretion through an accretion disk does not lead to a rapidly spinning SMBH, because 
either the accretion is intermittent (and residual coupling between the hole and surrounding 
gas can spin down the hole between the accretion episodes) or the accretion is not confined 
to a single plane over the long timescales (10® yrs for a 1 Mq yT~^ accretion rate onto a 
10® M0 hole) required to spin up the hole. Rather, SMBH's acquire a significant fraction 
of the maximal Kerr angular momentum only through the essentially instantaneous event 
of coalescence with a second black hole of similar mass. 

We consider mergers in a population of parent galaxies each of which contains a single, 
nuclear SMBH with small or zero spin. In a galaxy merger, the two nuclear SMBH's form a 
binary and eventually coalesce. The time scale for this coalescence is uncertain (Begelman, 
Blandford & Rees 1980), but we shall assume it is very much less than the Hubble time. 
Qualitatively, there are three possible outcomes of such merger. Two small (say ~ 10^ 
Mq) black holes of approximately equal mass will coalesce into a rapidly spinning, low 
mass hole. The merging of a large (say ~ 10® Mq) and a small hole will result in a slowly 
spinning, high mass hole. Lastly, the merging of two large holes will give a rapidly spinning, 
high mass hole. We propose that the radio galaxy and radio quasar phenomenon is a 
result of the last type of merger. The mass spectrum of SMBH's in the parent population 
may reasonably be expected to decline towards higher mass, so most mergers will involve 
relatively low mass holes and will not result in a high radio luminosity object. Because 
of the requirement that the holes are both massive and of comparable mass, the powerful 
radio galaxy/quasar phenomenon is intrinsically rare, and only results from a very small 
fraction of galaxy mergers. 
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In Section 2, we calculate the distributions of SMBH mass and jet luminosity for the 
merged population, while Section 3 discusses the inputs to the model in an astrophysical 
context. The results are compared with observations in Section 4 and discussed in Section 
5. Concluding remarks and some outstanding issues and problems are given in Section 6. 



We consider a parent population of galaxies each of which contains a slowly rotating 
or non-rotating SMBH of mass M in its nucleus. The number density of SMBH's in the 
parent population in the mass range M to M + dM is N{M)dM Mpc~^ in the mass range 
Ml to Mu- A fraction g of these galaxies is envisaged to merge with another galaxy in a 
given time and the associated black holes coalesce on a timescale -C the Hubble time. We 
neglect the loss of mass due to gravitational radiation and also assume that there is no 
spatial mass segregation of the parent population of SMBH's. 



We use the subscript m to represent parameters after merging. Two parents with 
masses and Mb combine to form a daughter with mass M^. The number density 
of merged SMBH's in the daughter population in the mass range Mm to Mm + dMm is 
Nm{Mm)dMm Mpc'^ with Nm{Mm) given by 



where No = Xvf" N{M)dM is the total number of SMBH's per Mpc^ of the parents, M2 
is the smaller of M„ and Mm — Mi, and Mi is the larger of Mi and Mm — M^. 

For a power-law distribution N{M) = KM~^ [K is a normalization constant) and 
substituting x = M/Mm, we have 



2. THE MODEL 



2. 1 Distribution of Mass after Merging 




(1) 




(2) 
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where X2 is the smaUer of My^/M^ and (M^ — Mi)/Mm, and xi is the larger of Mi/M^ 
and {Mm-Mu)/Mm- The daughter distribution in the mass range 2M/ < < M^ + Mi 
(i.e. a;i = Mi /Mm <^ 1/2, and a;2 = 1 — Mi /Mm — 1) does not show end effects of the 
mass range of the parents, and in this range we have 



Nm{Mm) = |^M-2^+ii?(l - 7, 1 - 7) for 7 < 1, 



(3) 



where B is the Beta function, and 



Nm{Mm)=gKM-^ 



for 7 > 1. 



(4) 



Thus the slope of the daughter mass function is identical to that of the parent if 7 > 1. 

The integral for Nm{Mm) may alternatively be expressed in terms of the mass ratio 
of the two merging black holes r = Ma/Mi, 

• 1 



Nm{Mm) = 



9M„ 



N 



ri(M^) 



MmV 

(1 + r) 



N 



(1+r) 



dr 



(5) 



where ri is the larger of Mi /{Mm - Mi) and {Mm - Mu)/Mu. We have taken r < 1 (M„ < 
Mfe) in order to avoid double-counting. For the power-law distribution N{M) — KM~'^ , 
equation (5) becomes 



1 -7 



{i + ry 



dr 



(1 + r) 



2 ■ 



(6) 



'ri(M^) 

Fig. 1 is a plot of lines of constant Mm and r in the Ma-Mi, plane. Each value of Mm 
describes a curve for which a specific range of r, corresponding to the range of integration 
in equation (5), is allowed. 

2.2 Distribution of Jet Luminosity after Merging 

The non-thermal emission is assumed to be powered by the spin energy of the merged 
SMBH. Theories of the physical processes by which power is extracted from the spinning 
hole are still in their infancy, so our calculation of the jet power must be considered prelim- 
inary and illustrative. It is believed that energy may be extracted from a spinning black 
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hole via the magnetic field by the Blandford-Znajek process (Blandford & Znajek 1977). 
The electromagnetic power available from this process depends on the hole's mass and spin, 
and the axial magnetic field generated by external currents. Under the assumption that 
the magnetic energy density near the black hole is proportional to the shear stress in the 
accretion disk, the resulting luminosity available for extraction is conjectured (Blandford 
1990) to be 

Lem = OLH^^LEdd- (7) 

where an is an unknown constant, LEdd is the Eddington luminosity, and A is the black 
hole spin parameter (see Blandford 1990 and references therein). In this formalism, the 
relationship between the "electromagnetic mass" (defined as = A^M^) and the 

jet luminosity (Lem = ctH-^^^^^-^em) is analogous to that between and the Edding- 
ton luminosity. We shall use equation (7) to obtain the jet luminosity, while recognising 
the idealised nature of the expression. 

We thus need to obtain the spin parameter A of the daughter black hole from the 
merging of the two parent black holes, which are assumed to possess no spin. The param- 
eter A is a measure of the spin per unit mass of the daughter and so, assuming a circular 
orbit for the pre-merged holes, A is completely determined by the ratio r of the masses 
of the parents. We may write A in terms of the specific angular momentum, a, and the 
gravitational radius, m = GMjc?^ of the SMBH (cf. Blandford 1990; in relativists units 
[G = c = 1], both a and m have units of length, so a/m is dimensionless) : 

a 

A= (8) 

The calculation of ^ from the coalescence of two black holes with arbitrary mass ratio is a 
difficult problem in general relativity that is currently unsolved. However, the functional 
forms for A and ^ can be constrained using physical arguments. For reasons stated in 
Appendix A, we have chosen the following parameterization for ^(r): 



m \ {1 + ry 



2\/3+ (3.95-2^3) 



(1+r) 
10 



(1+r) 



(9) 



where /i is the fraction of angular momentum at the last stable circular orbit which is lost 
to gravitational radiation during the plunge to coalescence for equal-mass parents. For 
small values of r, A(r) has the limiting form y/Sr. 

Sample plots of A(r) are shown in Fig. 2. We have experimented with other reasonable 
(see Appendix A) estimates of the form for A(r) when calculating the distribution of jet 
luminosities (see below), and note there are no significant effects on our results. We vary 
/i from 0.0 to 0.5 in the following calculations. 

The distribution of Mem after merging can be written 

NoJrr(X,Mem) A2(r)(l + r) X^{r){l + r) X^{r){l + ry 

The lower limit, ri, is determined by solving A^(ri)(l + ri) = Mgrn/My,. The upper limit, 
r2, is either 1 or the solution of A^(r2)(l + r2)/r2 = Mem/Mi, whichever is smaller. The 
changeover between these two values for r2 occurs at MeTn,crit = 2A^(1)M/. Fig. 3 shows 
lines of constant Mem and r in the Ma vs Mb plane. As in Fig. 1, each value of Mem 
describes a curve for which a specific range of r, corresponding to the range of integration 
in equation (10), is allowed. 

For a power-law parent mass function N{M) = KM~'^, equation (10) becomes 

NmiM^m) = ^M-i'^+' / r-^A^^-^(r)(l + O^^-^dr. (11) 

The integration in equation (11) can be performed analytically in limiting cases, and results 
are summarized in Appendix B. 

3. ESTIMATES OF IMPORTANT PARAMETERS 

In order to apply the model described in Section 2 and derive the properties of the 
merged ("daughter") population, we need to estimate a number of required parameters. 
We first describe how the parent mass function is obtained. Then, in order to guide the 
comparison of the model with observations (Section 4), we also estimate canonical values 
for g, an, and e (the efficiency of conversion of jet power to radio emission). 
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3.1 The Parent Mass Distribution 

This distribution may be estimated from the optical luminosity function (OLF) of 
AGN (e.g. Marshall 1985, 1986; Boyle, Shanks & Peterson 1988) by making assumptions 
about the bolometric correction (i.e. the ratio of total thermal luminosity to that in the 
B band), the Eddington ratio and the activity patterns of the AGN's over cosmological 
timescales. In particular, Cavaliere & Padovani (1988, 1989) have derived mass distri- 
butions under three assumptions concerning the activity patterns. At one extreme is a 
pattern (which they label 'C' for 'continuous') of literal luminosity evolution, with each 
object dimming continuously over the timescale (~ 3 Gyr) for evolution of the quasar 
population. High redshift quasars evolve to local Scyfert I's in a one-to-one relationship, 
and AGN activity is hosted by a privileged set of some 1% of bright galaxies. A second, 
intermediate pattern (called 'R' for 'recurrent' by Cavaliere & Padovani) is provided by 
discontinuous activity in a larger fraction of galaxies, consisting of short, recurrent episodes 
of activity. This behavior could be driven by, for example, episodes of accretion stimulated 
by interaction of the host galaxy with companion. In the third pattern (labelled 'E' for 
'Eddington'), which is at the opposite extreme to 'C', activity represents a single, short 
event in the lifetime of every galaxy. 

In pattern 'C', L/Le is estimated to be ~ 2 x 10~'^r]E,-i at the present epoch, where 
rjE,-i is the efficiency of mass - energy conversion, in units of 10~^. A small fraction of 
galaxies (i.e. those which exhibit activity) contain very massive (10^° -10^2 Mq) remnants. 
For the case of pattern 'R', Cavaliere & Padovani (1988) favor L/Le — 5 x 10~'^r]E,-i- The 
black hole masses for both currently active and inactive galaxies cover a range ~ 10^ — 10^ 
Mq , with currently inactive galaxies having a space density ~ 500 times that of currently 
active ones (see Fig. 1 of Cavaliere Sz Padovani 1988). For pattern 'E' in its simplest form, 
L/Le is taken to be ~ 1 and all galaxies contain SMBH with masses ranging between 
~ 10^ and 10^^ M©. To illustrate the application of our model, we shall use (Section 4) 
the Cavaliere Sz Padovani mass functions for cases 'C and 'R', derived using L/Le — 10~^ 
and L/Le — 5xl0~^, respectively. 
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3.2 The Merger Rate 



The parameter g represents the fraction of SMBH's in a unit volume of space that 
merge in the relevant time interval. Because the merged population of interest comprises 
those galaxies currently active as radio-loud AGN's, this time interval is the average lifetime 
of the extended radio sources, which is commonly estimated to be tr ~ 10^ years. The 
fraction of black holes which merge in this time depends on both the galaxy merger rate 
and the time required for the two SMBH's to coalesce. Both of these timescales are 
controversial. The merger rate is undoubtedly a function of both the environment and 
the mass ratio of the two galaxies. Schweizer (1982, 1990) argues that all ellipticals may 
have undergone at least one merger per Hubble time th- Merger rates of 1 per 10^ yrs 
may apply to close binary galaxies (Borne 1984) or groups (Navarro, Mosconi & Garcia 
Lambas 1987). Roos (1981) estimates a merger rate of 1 per 10^ yrs when the second 
galaxy has 1% of the mass of the first. Heckman, Carty & Bothun (1985) have shown 
that the average local galaxy density is larger (by at least a factor of 2 - 3) for radio-loud 
ellipticals and lenticular galaxies than for radio-quiet galaxies of the same morphological 
type. Therefore, it is likely that the somewhat higher merger rates (i.e. 1 in 10^ yrs) apply 
to radio galaxies. The time to coalescence, re, of the SMBH's has been variously claimed 
to be as short as ^ 10^ years (Ebisuzaki, Makino & Okumura 1991) or longer than the 
Hubble time (e.g. Valtaoja, Valtonen & Byrd 1989). In our model, we obviously require 
Tc « Th, and probably tc < 10^ yrs (so that the radio galaxy is recognisable as a recent 
merger). Writing tm = 10^ yrs for the average time between mergers, we take as our 
canonical value g ~ tr/tm — 10^/10^ ~ 0.1. 

For tr ~ O.Itm, it is unlikely that the merged galaxy will suffer a second merger 
during its lifetime as a radio source. Therefore, we do not need to consider the effect of a 
second black hole merger event on the mass and spin of the merged hole during the lifetime 
of most radio sources. 

3.3 The Electromagnetic Extraction Efficiency 
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The parameter an is a measure of the efficiency of conversion of the spin energy of 
the daughter SMBH to jet power. We shall adopt an = 1 as a canonical value. 

3.4 The Efficiency of Conversion of Jet Power to Radio Emission 

Here the typically quoted numbers of e = Lradio/Lem — 10~^ — 10~^ may be considered 
canonical (see Leahy 1991, who also discusses the limitations of this form of parameteri- 
sation of the radio source energetics). 

The effects of changes in these parameters on the properties of the OLF and radio 
luminosity function (RLF) of the merged (radio-loud) population are easily gauged. In- 
creases in N(M) (i.e. a greater density of parent holes of a given mass) and g increase 
the daughter number density at a given mass and luminosity. Increasing M increases the 
optical and radio luminosities, while altering an and e affects only the radio luminosities. 
In fact, the radio luminosities depend on the product an^, not on an and e separately. 

4. COMPARISON WITH OBSERVATIONS 

In comparing our model with the OLF and RLF of active galaxies, we have proceeded 
as follows. 

(a) A 'parent' population is first adopted. In the case of 'C type activity, only the currently 
active galaxies arc considered to contain SMBH's, so the appropriate parent is the radio- 
quiet population. Because the space density of radio-louds is much less than that of 
radio-quiets, we can, to a good approximation, use the total (radio-quiet plus radio-loud) 
population as the parent. For 'R' type activity, the whole galaxy population, which is 
dominated by currently inactive galaxies (see Fig. 1 of Cavaliere & Padovani 1988), should 
be taken as the parent. 

(b) To derive the parent mass distribution, we first fitted the B band OLF ^tot{MB) of 
the total population of AGN (Padovani 1993; Cheng et al. 1985) with a two power law 
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form (cf. Marshall 1986; Boyle, Shanks & Peterson 1988). The function chosen was: 



L0, which corresponds to = -23.65. The data and the fitted function are plotted in 
Fig. 4. By adopting a value, b, for the fraction of thermal emission which is radiated in the 
B band {b ~ 0.03; Blandford 1990), the total thermal luminosity corresponding to a given 
Mb rnay be found. Given the assumed Eddington ratios (L/Le) adopted by Cavaliere & 
Padovani (1988) and summarised above, the distribution of black hole masses in currently 
active galaxies is then calculated. For activity pattern 'C, this distribution represents the 
total mass distribution of SMBH's in galaxies, and is thus assumed to represent N{M). 
For activity pattern 'R', we must multiply the mass function by the ratio of the total 
number density of SMBH's (assumed to be dominated by currently inactive galaxies) to 
the number density of currently active SMBH's. We have used the ratio ~ 500 advocated 
by Cavaliere & Padovani (1988, 1989) to obtain N{M) for pattern 'R'. An important 
issue is the relationship between the activity patterns of the radio-quiet and radio-loud 
populations. For 'C type activity, we assume that a radio source turns on shortly after 
the SMBH's have merged and continues as a radio-loud AGN for a single period of order 
10^ yrs (cf. Section 3.2). For 'R' type activity, the radio-quiet and radio-loud patterns 
are assumed to be similar: only 1 in 500 of the merged holes is radio active at any given 
time and, when active, the lifetime of the activity is the canonical 10^ yrs. Physically, the 
inactive times would correspond to periods without gaseous accretion; such accretion is 
necessary both for fuelling the thermal emission and for providing the magnetic field used 
to tap the spin energy of the hole. 

(c) The distribution of merged mass, N^iMm), was then calculated using equation (5). 
In fact, the ratio Nm{Mm)/N{M) = g when Mm = M (equation 4). If we assume that 
the values of h and L/Le are the same for the radio-loud and radio-quiet objects, then 




(12) 



where TV* = 5 x 10 ^ Mpc ^. The luminosity function changes slope at L*^ = 5.8 x 10 
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Nri{MB)/Nrq{MB) = 9, independent of b and L/Le- Here Nri{MB) is the OLF of the 
radio-louds and Nrq{MB) that of the radio-quiets. By fitting our calculated Nj.i{Mb) to 
Padovani's measurements, we obtain ^ = 0.1 (see Fig. 4), although the slightly higher 
value ol g = 0.3 may apply at high optical luminosities {Mb < —24). For a radio galaxy 
lifetime of ~ 10^ yrs, the implied rate of merging is 1 in 10^ yrs, in excellent agreement 
with that obtained from other considerations in Section 3.2. 

(d) Given the parent mass distribution N{M) (see (b) above), the value of 51 just obtained, 
and the assumed forms of A(r) (equations 8 and 9), we calculate Nm{Mem) using equation 
(10). N(Lern) is then found (using Lgm = an ^^^^^ ^ -^em ; we have taken an = 1), with 
Lem assumed to represent the jet luminosity. Adopting a value for the efiiciency, e, of 
conversion of jet power to radio emission and a form for the synchrotron spectrum (a 
power law S oc with a = 0.75 between 10 MHz and 100 GHz), we finally obtain 
the RLF of radio- loud AGN in terms of P20, the monochromatic radio power at 20 cm. 
The predicted curves are plotted along with the observationally determined RLF's of radio 
galaxies and radio-loud quasars in Fig. 5. 

5. DISCUSSION 

We have found that our model can account for the present day relative numbers of 
radio-loud and radio-quiet AGN's if the merger rate is 1 in 10^ yrs. However, if such a high 
rate were sustained over a Hubble time and is applicable everywhere, the parents (spirals) 
would have long since disappeared. It is known that the galaxies hosting radio-quiet AGN's 
tend to be in low galaxy density environments whereas those hosting radio-loud objects 
are found preferentially in high density environments (e.g. Smith & Heckman 1990). A 
more sophisticated calculation would allow for a range of surrounding galaxy densities and 
hence merger rates; radio-quiet (unmerged) objects would naturally tend to be found in 
low density and radio-loud (merged) in high density regions. Elliptical-spiral mergers are 
likely to be important in the denser regions. 

The good fit to the RLF shown in Fig. 5 is obtained through mergers involving the 
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low mass (M^ > -24) portion of the parent mass function (Fig. 4). In this region, the 
mass function of the parents is N(M) oc M~^-^, which imphes a RLF A^rn(-p2o) cx: -P2^^'^ 
(equation B4), in good agreement with the RLF observed above the break near logP2o = 25 
WHz~^. Mergers involving the high mass part of the parent mass function, which follows 
N{M) cx M~^-^, are predicted (equation B4) to give Nm{P2o) oc -P20^' ^ much steeper 
form than observed in any part of the RLF. This steep predicted form at high masses is 
intrinsic to the conjecture (equation 7) used to calculate the radio luminosities. The formal 
implication is that black hole masses in FRII radio galaxies should be similar to those in 
radio-quiet objects with Mb > —24, rather than those with Mb < —24. 

If the parent mass distribution, N{M), has a sharp cut-off at masses below those 
probed by the OLF (Fig. 4), the predicted RLF function for logP2o < 10^^ WHz'^ 
(corresponding to Mem < Mern,crit — 2A^(1)M2) has the form Nm{P2o) oc P^o^ (equation 
B5). This last form does not fit the data near or below the break in the RLF, suggesting 
there is a low luminosity turnover, rather than a sharp cut-off, in the OLF at Mb > —18. 
Alternatively, it would be possible to make an ad hoc modification of the form of equation 
(7) to allow the break in the OLF (i.e. the parent mass function) to correspond to the 
break in the RLF. 

In the context of the current model, cosmological evolution of the radio luminosities 
would most naturally be ascribed to a higher black hole merger rate at earlier epochs. 
Such would provide density evolution, but radio luminosity evolution could also occur if 
an and/or e are epoch dependent. 

6. CONCLUDING REMARKS AND FUTURE WORK 

In this paper, we have assumed that radio-loud AGN have rapidly spinning, super- 
massive black holes (SMBH's) in their nuclei and that the radio jets are powered by the 
spin energy of the holes. Contrary to the conventional view that SMBH's are spun up 
via gaseous accretion through a disk, we argue that such holes acquire a large spin only 
through a merger with another SMBH. Presumably low mass SMBH's are more common 
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in galactic nuclei than high mass ones, so most mergers will involve relatively low mass 
SMBH's. Let us also suppose that the mass of the parent galaxy is correlated with the mass 
of its nuclear SMBH. Bearing in mind that only mergers between galaxies of approximately 
equal mass are likely to give elliptical-like remnants, we may qualitatively distinguish three 
types of mergers, depending on the masses of the two SMBH's (and hence parent galaxies) 
involved. 

(i) Low mass - low mass coalescences will lead to a rapidly spinning, low mass hole in 
an elliptical-like galaxy. Such objects are probably those ellipticals with low luminosity, 
pc-scale radio cores, which are common in early-type galaxies (e.g. Sadler et al. 1994). 

(ii) Low mass - high mass coalescences result in a slowly spinning, high mass hole. The 
merged galaxy will retain the morphological type of the more massive galaxy and such 
galaxies will form part of the radio-quiet population. However, these objects are expected 
to comprise only a small fraction of the radio-quiet population, which is dominated by the 
unmerged galaxies (parents). 

(iii) High mass - high mass coalescences form a rapidly spinning, high mass hole. We 
associate these rare events with the formation of a powerful radio source. 

Our idealised calculations show that both the optical and radio luminosity functions 
of powerful radio galaxies can be reproduced by a simple model. Merger rates of 1 per 
10^ yrs are required in the regions where FRH radio galaxies form. The product of the 
parameters an (defined by an — Lem/ ^'^ir)LEdd, where Lem is the jet luminosity, A the 
black hole spin parameter, and LEdd is the Eddington luminosity) and e (e = Lradio/Lem) 
is constrained to be ^ 5xl0~^ - lO"'^, depending on which pattern of evolution of AGN 
activity is chosen (Section 4) . 

In summary, our model accounts naturally for the following observations: 

(1) No spiral galaxies host radio- loud AGN. 

(2) Powerful radio galaxies and quasars are the products of a recent merger. 

(3) Only a tiny fraction of galaxy mergers produce a radio galaxy (despite the belief that 
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SMBH's are common in galactic nuclei). 

(4) Radio- loud AGN are much rarer than radio-quiet AGN. 

We conclude by listing a number of issues and problems which require further work. 

(a) Measurements of the galaxy density and merger rate in the environment of powerful 
radio sources are needed. Further observational work on FRII radio galaxies is particularly 
desirable. It is, of course, possible that some radio galaxies are the result of the merger 
of an isolated binary galaxy. In that case, the density of galaxies around the radio galaxy 
after the merger would be a poor indicator of the relevant merger rate. 

(b) Further calculations of the timescale for coalescence of two nuclear SMBH's in a galaxy 
merger due to stellar and gas dynamical processes are important. 

(c) The controversial (e.g. Punsly & Coroniti 1990) topic of energy extraction from spin- 
ning SMBH's deserves more study, in order to provide a firmer estimate of the jet power. 

We thank R. D. Blandford, J. P. Henry, A. lUarionov, T. Jacobson, L. Kidder, M. 
Livio, E. Poisson, and B. C. Whitmore for valuable discussions. This work was supported 
in part by NASA grants NAGW-2689 and NAGW-3268. 



19 



REFERENCES 



Antonucci, R. R. J. 1993, ARAA, 31, 473 

Auriemma, C, Perola, G. C, Ekers, R., Fanti, R., Lari, C, Jaffe, W. J. & Ulrich, M. H. 

1977, A&A, 57, 41 
Baum, S. A. & Heckman, T. M. 1989, ApJ, 336, 702 
Begelman, M. C, Blandford, R. D. & Rees, M. J. 1980, Nature, 287, 307 
Blandford, R. D. 1990, In Active Galactic Nuclei, ed T. J.-L. Courvoisier & M. Mayor, 

Saas-Fee Advanced Course 20 (Berlin: Springer- Verlag) , 161 
Blandford, R. D. & Znajek, R. L. 1977, MNRAS, 179, 433 
Borne, K. D. 1984, ApJ, 287, 503 

Boyle, B. J., Shanks, T. & Peterson, B. A. 1988, MNRAS, 235, 935 
de Bruyn, A. G. & Wilson, A. S. 1978, A&A, 64, 433 
Cavaliere, A. & Padovani, P. 1988, ApJ, 333, L33 
Cavaliere, A. & Padovani, P. 1989, ApJ, 340, L5 

Cheng, F., Danese, L., dc Zotti, G. & Franceschini, A. 1985, MNRAS, 212, 857 
Ebisuzaki, T., Makino, J. & Okumura, S. K. 1991, Nature, 354, 212 
Edelson, R. A. & Malkan, M. A. 1986, ApJ, 308, 59 

Elvis, M., Green, R. F., Bechtold, J., Schmidt, M., Neugebauer, G. Soifer, B. T., Matthews, 

K. & Fabbiano, G. 1986, ApJ, 310, 291 
Fanti, R. & Perola, G. C. 1976, In Radio Astronomy & Cosmology, lAU Symposium Nr. 

74, ed D. L. Jauncey pl71 (Dordrecht: Reidel) 
Heckman, T. M., Carty, T. J. & Bothun, G. D. 1985, ApJ, 288, 122 
Heckman, T. M. et al. 1986, ApJ, 311, 526 
Hutchings, J. B. 1987, ApJ, 320, 122 

Hutchings, J. B., Janson, T. & Neff, S. G. 1989, ApJ, 342, 660 
Hutchings, J. B. & Neff, S. G. 1990, AJ, 99, 1715 

20 



Kellermann, K. I., Sramek, R. A., Schmidt, M., Shaffer, D. B., & Green, R. 1989, AJ, 98, 
1195 

Kidder, L. E., Wiii, C. M., & Wiseman, A. G. 1993, Physicai Review D, 47, 8, 3281 
Kormendy, J. & Djorgovski, S. 1989, ARAA, 27, 235 

Leahy, J. P. 1991, In Beams and Jets in Astrophysics, Ed. Phiiip A. Hughes, p 100 

(Cambridge University Press) 
Maikan, M. A. 1984, ApJ, 287, 555 
Marshaii, H. L. 1985, ApJ, 299, 109 

Marshaii, H. L. 1986, In Structure and Evoiution of Active Galactic Nuclei, Eds G. Giuricin 

et al., p 627 (Reidel: Dordrecht) 
Meurs, E. J. A., & Wilson, A. S. 1984, A&A, 136, 206 
Miller, L., Peacock, J. A., & Mead, A. R. G. 1990, MNRAS, 244, 207 
Navarro, J. F., Mosconi, M. B., & Garcia Lambas, D. 1987, MNRAS, 228, 501 
Neugebauer, G., Miley, G. K. Soifer, B. T. & Clegg, P. E. 1986, ApJ, 308, 815 
Padovani, P. 1993, MNRAS, 263, 461 

Peacock, J. A., Miller, L., & Longair, M. S. 1986, MNRAS, 218, 265 
Punsly, B. & Coroniti, F. V. 1990, ApJ, 350, 518 

Rawlings, S. & Saunders, R. 1991, Nature, 349, 138 

Rawlings, S. 1994, In The First Stromlo Symposium: The Physics of Active Galaxies, ASP 
Conference Series Vol. 54, Eds G. V. Bicknell, M. A. Dopita, & P. J. Quinn, p 253 
(Astronomical Society of the Pacific: San Francisco) 

Rees, M. J. 1984, ARAA, 22, 471 

Roos, N. 1981, A&A 104, 218 

Sadler, E. M., Slee, O. B., Reynolds, J. E. & Ekers, R. D. 1994, In The First Stromlo 
Symposium: The Physics of Active Galaxies, ASP Conference Series Vol. 54, Eds 
G. V. Bicknell, M. A. Dopita, & P. J. Quinn, p 335 (Astronomical Society of the 
Pacific: San Francisco) 

21 



Sanders, D. B., Phinney, E. S., Neugebauer, G., Soifer, B. T. & Matthews, K. 1989, ApJ, 
347, 29 

Schmidt, M. 1970, ApJ, 162, 371 

Schweizer, F. 1982, ApJ, 252, 455 

Schweizer, F. 1990, Talk given at STScI Workshop, "The Galaxy Merger Rate" 
Smith, E. P., Heckman, T. M., Bothun, G. D., Romanishin, W. R., & Balick, B. 1986, 
ApJ, 306, 64 

Smith, E. P. & Heckman, T. M. 1990, ApJ, 348, 38 
Steidel, C. C. & Sargent, W. L. W. 1991, ApJ, 382, 433 

Strittmatter, P. A., HiU, P., Pauliny-Toth, I. I. K., Steppe, H., & Witzel, A. 1980, A&A, 
88, L12 

Toomre, A., & Toomre, J. 1972, ApJ, 178, 623 
Valtaoja, L., Valtonen, M. J., & Byrd, G. G. 1989, ApJ, 343, 47 
Veron-Cetty, M. P. & Woltjer, L. 1990, A&A, 236, 69 
Wall, J. v., Pearson, T. J. & Longair, M. S. 1980, MNRAS, 193, 683 
Woltjer, L. 1990, In Active Galactic Nuclei, ed T. J.-L. Courvoisier & M. Mayor, Saas-Fee 
Advanced Course 20, pi (Berlin: Springer- Verlag) 



22 



Figure Captions 



Figure 1. Contours of the merged mass Mm (solid) and the mass ratio r (dashed) in the 
Ma-Mb plane for the case M/ = 10^ M©, M„ = 10^ M0. In order to avoid double-counting, 
we have chosen the upper limit of integration in r (r < 1) from the condition Ma < Mf,. 
The dotted contour is = My^+Mi. When Mm < My^+Mi, the lower limit of integration 
in equations (5) and (6) is determined from Ma = Mi, i.e. ri = Mi/ [Mm — Mi). When 
Mm > Mu-\- Ml, the lower limit is determined from Mb = M„, i.e. ri = {Mm - Mu)/Mu. 

Figure 2. The spin parameter A as a function of the mass ratio r for /i = 0.0 - 0.5. The 
limiting form for small r (A(r) = \/3r) is shown as a dashed line. 

Figure 3. Contours of the merged electromagnetic mass M^m (solid and dotted) and the 
mass ratio r (dashed) in the Ma-Mf, plane. The solid lines represent constant M^m for fi 
= 0, while the dotted represent constant Mem for fi = 0.5. The lower limit of integration 
(ri) in equations (10) and (11) is determined from M^ = Mu- The upper limit (r2) is 
determined by M„ = M^ {r = 1) for M^m > Mem,crit = 2A2(l)Mi, and by M« = Mi for 

Figure 4. Optical Luminosity Functions of Active Galaxies. Filled squares represent 
all quasars in the Padovani (1993) sample, while filled triangles correspond to Seyfert 
galaxy nuclei (Cheng et al. 1985). Open squares correspond to the radio-loud quasars 
in Padovani's sample. The parent mass function was computed by fitting the luminosity 
function for the whole sample (upper line) with a two power-law form (see Section 4). 
The merger rate oi g = 10~^ was chosen so that the predicted luminosity function for 
the daughter (merged) population (lower curve) agrees with the data for the radio-loud 
quasars (Section 4). 

Figure 5. Radio Luminosity Functions of Active Galaxies for /i = 0.0 - 0.5. The sym- 
bols represent the observationally inferred radio luminosity functions for radio galaxies 
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(Auriemma et al. 1977 - squares), radio galaxies plus radio-loud quasars (Wall, Pearson 
& Longair 1980 - triangles) and radio- loud quasars (Fanti & Perola 1976 - circles). The 
dashed lines represent extensions of the parent mass function beyond the range in Mb of 
the Cheng et al. (1985) and Padovani (1993) samples (sec Figure 4). For the 'C type 
models, the plotted curves correspond to e = 10~^, while for the 'R' type, e = 5xl0~^. 

Figure Al. The spin ^ as a function of the mass ratio r for fi = 0.0 - 0.5. The limiting 
form for small r is shown as a dashed line. 

Figure A2. The derivative of the spin parameter A with respect to the mass ratio r as a 
function of r for /i = 0.0 - 0.5. 
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Appendix A: The Spin Parameter A 

The spin parameter A and ^ are related through the transformation 



The problem of finding the resulting spin ^ as a function of r, the ratio of the masses of 
the two parent black holes, involves general relativistic computer codes and a solution is 
not yet available. A recent study by Kidder, Will & Wiseman (1993) evolves the binary 
system down to the innermost stable circular orbit (ISCO) for arbitrary mass ratio. They 
find (see their Fig. 3) that the angular momentum per unit reduced mass (J) at ISCO 
is approximately linear in the parameter rj (the ratio of the reduced mass to the total 
mass), rising from ^ = 2^/3 for the test-mass case (r— >0, r/— >0) to ^ = 3.95 in the case 
of equal masses (r = 1, ry = 0.25). The parameter ij can be written in terms of the mass 
ratio as ?7 = r/(l + r)^. In order to bridge the gap between ^ at ISCO and ^ after 
coalescence, we introduce a function / (the fraction of angular momentum at ISCO that is 
lost in gravitational radiation during the plunge from ISCO to coalescence), which is also 
a function of the mass ratio r. Then we can write 

^r) = (l-/(r))f^^ ry (^2) 



m \my jsco 



(1 - m) 



2^3 +(3.95 -2^3)- 

(1 + r) 



4r 

~2 



(l + r)2 

We may place restrictions on /, ^ or A in order to constrain the functional form A(r). 
We must have 

< — (r) < 1 and < A(r) < 1 for < r < 1. (i) 
m 

It is obvious that we must also have 

< /(r) < 1 for < r < 1. (ii) 



With this restriction, the requirement that ^ < 1 always holds true. 
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In the test- mass limit, we assume that the loss of angular momentum during the 
plunge is negligible (/(O) = 0), which yields a limiting form for A(r): 

X{r) = ^ — (r) = VSr as r — > 0. (in) 

One expects the spin to increase monotonically with the ratio of the masses of the 
parent black holes, which implies an additional requirement of 

> 0, or, equivalently, ^>0 0<r<l. (iv) 
dr dr 

Various functional forms for /(r) were tried. A simple linear form /(r) = a + br 
(a =0, b < 1) fails criterion (iv). Other forms {e.g. polynomial, polynomial/exponential 
combination) can be constructed to satisfy (i — iv), but such forms require either a large 
number of parameters (in which case interpreting variations in a single parameter is not 
straightforward) or a complex functional form. A linear function in r] [f (r) — a + br]; a = 
0, 6 < 4) does satisfy (i — iv) if b <2.1, so we have used the single-parameter form 

in our calculations (section 2.2). Here fi = b/4: = /(I) is the fraction of angular momentum 
at ISCO lost to gravitational radiation during the plunge in an equal-mass system. The 
resulting spin ^ and spin parameter A are shown in Figs. Al and 2, respectively, for values 
of /i = 0.0 - 0.5. For all values of /i, ^ = \/3 at r = and > for aU r (see Fig. A2). 
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Appendix B: Analytic Forms for the Distribution of Mgm 

Equation (11) can be simplified in limiting cases. In the following, we have 
A(r) = \/3r for all r: 

7 < 1 and 3Mf /M„ < Mem < Mu : 



„^2327-3 I 3 -(-y+i) 



7=1 and Mgm ,crit 

< Mem < 



7 = 1 and ZMf/Mu < Mg^ < Me^,crit 



Nm{Mem) = Me„^ ^^M^ V 3 ^" 



7 > 1 and Mem ,crit 

< Mem < Mu : 



N (M )- ^^^^-^M "^'^+^ 



'em 



7 > 1 and 3Mf/Mu < Mg^ < Mgm.crit 



AT CM W ^ f_!_)37-3j^7-2 
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